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I.  SUMMARY 

Operation  of  the  monopropellant  rocket  motor  has  been  performed  with 

artificial  flow  modulation  at  frequencies  up  to  200  cycles  par  second  at 

several  chamber  pressures,  and  some  satisfactory  measurements  of  tne  pres* 

* 

sure  lag  between  injector  and  chamber  have  been  made.  Eventually,  this 
data  should  demon strata  the  suitability  of  the  tocanique  employed  to  make 
measurements  and  at  the  same  time  provide  a rough  check  on  the  fundamental 
precepts  of  Crocco's  theoretical  analyses  as  presented  in  earlier  reports0 
Improvements  in  the  accuracy  of  the  measurements  and  both  theoretical  and 
experimental  corrections  to  the  instantaneous  flow  rate  determinations 
must  be  completed  before  the  data  can  be  fully  analyzed# 

Several  methods  of  measuring  the  pressure  lag  have  been  evolved 

utilizing  DC  amplifiers,  constant-phase  AC  amplifiers,  and  direct  recording 

of  pressure  traces  on  a four-channel  oscilloscope.  Amplitude  calibrations 

have  not  yet  achieved  the  desired  degree  of  accuracy,  but  it  is  expected 

to  be  only  a matter  oi'  time  until  auditionai  refinement  of  circuit  components 

% 

will  permit  obtaining  the  required  precision# 

The  Mittelmann  electromagnetic  flowmeter  proved  unsatisfactory  due  to 
equipment  malfunctions,  despite  the  many  hours  devoted  to  its  operation  and 
repair#  It  will  be  returned  to  the  vendor  for  additional  repairs#  Tne  ii 
mass  flowmeter  has  been  received,  and  preliminary  checks  indicate  that  the 
leakage  ha3  been  reduced  to  the  point  at  which  it  can  be  used  in  the  flow 
calibration  rig#  kul3ing  flow  tests  utilizing  this  motor  and  the  hot  wire  in 
conjunction  with  pressure  pickups  will  be  made  after  concluding  the  current 
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series  of  monopropeliant  rocket  motor  runso 

Detailed  information  on  the  performance  of  the  improved  Li-iiu 

v 

pressure  pickup  under  high-temperature  conditions  ha3  been  obtained  at  our 
request  by  the  NACA  Lewis  Flight  Propulsion  Laboratory*  heat  transfer  rates 
up  to  about  7o6  BTU  per  square  inch  per  second  (refered  to  the  pickup 
diaphragm  area)  have  been  experience^ aider  "screaming"  conditions  when  the 
pickup  was  flush-mounted  in  an  acid-gasoline  motor,  and  the  pickup  drift 
was  shown  to  be  anall  at  temperatures  up  to  $00°F«  *11  available  data  on 

this  pickup  design  are  now  being  collected,  and  will  be  presented  at  the 
Rocket  Combustion  Instability  Symposium  scheduled  for  28  October  1953* 

Construction  of  the  bipropellant  test  stand  is  about  one-third  completed, 
and  is  proceeding  according  to  schedule.  Details  of  the  system  are  presented 
later  in  the  report. 
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II  I.Ek'JJCTIJ.i 
~'io  object 

duAer  Contract  Joas  62-713-c  was  undortak on  as  a part  oi'  tno  jet 
propulsion  research  program  of  the  department  of  aeronautical  engineer- 
ing at  drinceton  to  "conduct  an  investigation  of  the  ;jonerai  proolen  of 
combustion  instability  in  liquid  propellant  rocket  engines 0 This  pro- 
gram shall  consist  of  theoretical  anal;  sea  and  experimental  verification 
of  tneory.  The  ultimate  objective  shall  be  the  collection  of  sufficient 
data  that  shall  permit  the  rocket  engine  designer  to  produce  power  plants 
'..■hich  are  relatively  free  of  the  phenomena  of  instability*  Interest  shall 
center  in  that  form  of  unstable  operation  which  i3  characterized  by  high 
frequency  vibrations  and  is  commonly  known  as  •screaming • y 0 

B.  HISTuHY 


intercut  at  Princeton  in  the  problem  of  combustion  instability  in 
liquid  proncilant  rocket  motors  was  given  impetus  by  a bureau  of  Aero- 
nautics symposium  held  at  the  ‘'‘aval  Resource  laboratory  on  tlio  ?th  and 
6th  of  nc cumber  This  interest  resulted  in  theoretical  analyses 

by  iTofessors  ii*  burjnerfieid  and  Lb  irocco  of  this  Center. 

Professor  Sunro  erf  i eld's  work, "Theory  of  Unstable  Comoustion  in 
liquid  Propellant  Rocket  Systems"  (JARS*  Sopt„  1R51),  considers  the 
effects  of  tooth  inertia  in  the  liquid  propellant  feed  lines  and  combus- 
tion chamber  capacitance  vith  a constant  combustion  time  lag,  :,nd 
applies  to  the  case  of  low  (up  to  about  200  cycles  per  second)  fre- 
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II  LJTRODUCTIuN  (Cont ' d. ) 

1 i»  History  (Cont'd.) 

quency  oscillations  sometimes  called  "chugging". 

Professor  Crocco  advanced  the  concept  of  the  pressure  dependence 
of  the  time  lag  in  mid-1951j  his  paper,  "Aspects  of  Combustion  Stability 
in  liquid  Propellant  Rocket  Motors"  (JARA,  Nov.  1951  and  Jan.-Feb.l952), 
presents  the  fundamentals  resulting  from  this  concept,  and  analyses  the 
cases  of  low  frequency  instability  with  monopropellants,  low  frequency 
instability  with  bipropeliants  and  high  frequency  instability,  with  com- 
bustion concentrated  at  the  end  of  the  combustion  chamber. 

Desiring  to  submit  the  concept  of  a pressure  dependent  time  lag  to 
experimental  test  a preliminary  proposal  was  made  by  the  University  to 
the  bureau  of  Aeronautics  in  ti»  summer  of  l?5l  and,  following  a formal 
request,  a revised  proposal  was  submitted  which  resulted  in  Contract 
Noas  52-713*c» 

Analytical  studies  with  concentrated  and  distributed  combustion  had 
been  carried  on  in  the  meantime  under  Professor  Crocco' s direction  and 
within  the  sponsorship  of  the  Uuggenheim  Jet  Propulsion  Center  by  S.  I. 
Cheng  and  were  issued  as  ids  Ph.D.  tuesis,  "intrinsic  High  Frequency 
Combustion  Instability  in  a liquid  Propellant  Rocket  Motor",  dated  April 
1952. 


Time  was  devoted,  in  anticipation  of  the  contract,  during  the  first 
third  of  1952,  to  constructing  facilities,  securing  personnel,  and  plan* 
ning  the  experimental  approach. 
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II  INTRODUCTION  (Cont'd.) 

B.  History  (Cont'd.) 

During  the  first  throe  month  period  of  the  contract  year,  personnel 
and  facilities  at  the  new  James  T’orrestal  Research  Center  were  assigned, 
and  the  initial  phases  of  the  experimental  program  were  planned  in  some 

detail* 

A constant  rate  monopropellant  feed  system  was  designed  and  pre- 
liminary designs  of  tne  ethylene  oxide  rocket  motor  and  the  instrumenta- 
tion systems  were  worked  out.  Special  features  of  the  projected  systems 
included  a pulsing  unit  to  cause  oscillations  in  propellant  flow  rate,  a 
water-cooled  strain-gauge  pressure  pickup  designed  for  flush  mounting  in 
the  rocket  chamber,  and  several  possible  methods  for  dynamic  measurement 
of  an  oscillating  propellant  flow  rate* 

Searches  were  made  of  the  literature  for  sources  of  information 
on  combustion  instability  and  ethylene  oxide,  and  visits  to  a number 
of  activities  working  on  iiquid  propellant  rocket  combustion  instability 
problems  were  made  for  purposes  of  familiarisation  witn  equipment  and 

results* 

The  basic  precepts  of  Crocco's  theory  for  combustion  instability 
were  reviewed,  and  detailed  analyst*  made  for  specific  patterns  of  com- 
bustion distribution* 

Operational  tests  and  calibration  of  the  i Vince ton- KIT  pressure 
pickup  proved  the  value  of  the  design ^ although  failure  of  the  pickup 
under  "screaming”  rocket  conditions  showed  the  necessity  for  modification 
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II  INANIMATION  (Contra) 

13 a History  (Cont*d«) 

of  the  cooling  system# 

Construction  of  the  monopropellant  test  stand  and  rocket  motor  was 
completeda  Modifications  were  made  to  the  Princeton-MIT  pressure  pickup 
to  provide  for  higher  permissible  heat-transfer  rates  in  order  that  it  be 
satisfactory  for  use  under  "screaming"  conditions  in  a bipropellant  rocket 
motor a Construction  and  preliminary  testing  of  the  hot-wire  flow  phas^ 
meter  and  its  associated  equipment  were  completeda 

A new  contract,  HKtos  53-817-c,  dated  1 March  1?53>  was  granted  by  the 
Bureau  of  Aeronautics  to  continue  the  program  original ly  started  under 
NOas  52-713-c.  Operation  of  the  monopropellant  rocket  motor  was  begun 
under  this  new  contract,  and  shakedown  operations  were  completed*  It  was 
found  that  decomposition  of  ethylene  oxide  could  not  be  attained  with  the 
original  motor  design  despite  many  configuration  changes,  and  it  was  deci- 
ded to  avoid  a long  and  costly  development  program  by  operating  the  "mono- 
propellant"  motor  with  small  amounts  of  gaseous  oxygen.  The  required 
limits  of  oxygen  flow  rate  were  determined  at  several  chamber  pressures, 
and  it  was  demonstrated  that  the  oxygen  would  probably  have  a negligible 
effect  on  performance  when  compared  to  the  effect  of  ethylene  oxide  flow 
rate  modulation*  Preliminary  tests  with  flow  rate  modulation  up  to 
100  ops  were  performed  for  the  purposes  of  system  checkout,  using  interim 
AC  amplifiers  in  lieu  of  the  necessary  DC  instruments* 

The  time  constant  of  the  hot-wire  liquid  flow  phasemeter  was  found  to 
be  0;15  milliseconds  and  preparations  for  instantaneous  flow  calibrations 
were  made*  A test  rig  was  constructed  for  this  purpose* 


ft*#  U.  Aero#  hngg#  Ri>t«  No,  2l6e 


10 


II  INFORMATION  (Cont«d.) 

B«  History  (Cont’d, ) 

A bipropellant  rocket  system  using  liquid  oxygen  and  IDO,,  ethyl 
alchbol  was  designed  on  the  basis  of  monopropellant  operational  experience 
incorporating  an  adjustable-phase  flow  modulating  unit  operating  on  both 
propellants.  Injector  design  was  based  on  a configuration  used  extensive- 
ly by  Reaction  Motors,  Inc, 

Subsequent  efforts  to  date  are  presented  in  detail  in  this  report. 
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III  APPARATUS 

A.  Monopropellant  Rocket  Motor 

Operation  of  all  components  of  the  motor  proper  continued  satisfact- 
orily during  this  report  period)  using  the  oxygen  flow  rates  indicated  in 
ths  previous  report. 

— * 

The  slowdown  of  the  flow  modulating  unit  experienced  on  earlier  runs 
was  traced  to  the  extremely  high  bearing  friction  which  occurred  when 
fluid  pressure  was  applied  to  the  piston*  rather  than  to  the  use  of  an 
underpowered  drive  motor  as  was  first  suspected.  Addition  of  a pressur- 
ised lubrication  system  to  the  laain  connecting-rod  sleeve  bearing  did 
not  correct  the  situation*  so  a new  split  shaft  was  manufactured  which 
permitted  the  use  of  a roller  bearing  in  this  application.  Furthermore* 
the  two 80- pound  flywheels  were  removed  from  the  unit.  This  was  originally 
done  in  order  to  reduce  the  excessive  time  required  for  the  shaft  to  attain 
constant  speed  (over  one  minute  at  6,000  RPM*  and  well  over  three  minutes 
at  12,000  RPM).  No  visible  effect  on  the  resulting  pressure-wave  shape 
was  observed  with  tils  configuration*  and  since  the  reduction  in  speed- up 
time  (10  seconds  at  12*000  RPM)  resulted  in  appreciable  reduction  in  mar 
of  critical  parts*  the  "flywheel* less"  unit  was  used  to  obtain  all  subse- 
quent test  data. 

B.  Bipropellant  Rocket  Motor 

ever all  design  of  the  bipropellant  feed  system  was  based  partly  on 
experience  obtained  with  the  nonopropellant  system  and  partly  on  general 
experience  gained  in  the  rocket  industry.  Although  the  first  propellant 
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III  APPARATUS  (Cont'd.) 

B.  Bipropellant  Rocket  Kotor 

combination  to  be  used  will  be  liquid  oxygen-100;j  ethyl  alcohol,  the  system 
was  designed  with  a view  toward  future  use  of  any  of  the  propellants  now 
in  common  usage  throughout  the  country,  with  the  possible  exception  of 
liquid  fluorine. 

Figure  1 shows  the  overall  schematic  layout  of  the  feed  system*,  The 
fuel  tank  is  a surplus  stainless  steel  vessel  of  approximately  20  gallon 
capacity,  hydrostatically  tested  to  2ii>00  psi0  A surplus  UO-gallon  Konel 
tank  with  a vacuum-in sulatod  jacket  was  also  obtained  and  hydrostatically 
tested  to  1500  psi.  for  use  with  liquid  oxygen  at  1000  psi  maximum  pressure* 
This  latter  tank  will  be  suitable  for  the  first  few  months  of  operation, but 
testing  at  the  higher  chamber  pressures  will  require  a vessel  of  greater 
strength.  Hie  vent-relief  valves  are  of  the  same  design  as  the  one  used  so 
satisfactorily  on  the  uonoprop ellant  tank,  but  are  constructed  of  stainless 
steel;  with  the  liquid  oxygen  valve  protected  from  malfunction  due  to  frost- 
ing. A pressure- actuated  emergency  propellant  valve  from  too  Riil  UOO-pound 
thrust  Lark  motor  has  been  modified  for  use  a:;  the  main  propellant  valve? 
Cavitating  venturis  of  three  different  flow  ratings  will  be  used  for  the 
tliree  chamber  pressure  levels  anticipated,  and  Potter  flowmeters  will  again 
be  used  to  reoord  stoady-stato  flow  rates* 

A flow  modulating  unit  for  operation  at  frequencies  up  to  200  cycles 
per  second  has  been  constructed  for  use  in  this  system.  The  unit  consists 
of  two  easily  replaceable  piston-cylinder  assemblies?  with  the  pistons 
driven  from  a coman  shaft  (see  Figure  2).  In  order  to  avoid  variation  in 
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III  APPARATUS  (Cont'd. ) 

B.  Bipropellant  Socket  iiotor  (Cont'd) 

propellant  mixture  ratio,  the  crank  angle  between  the  two  piston  strokes 
is  continuously  adjustable  from  zero  to  360°  so  that  simultaneous  modula- 
tion of  the  two  propellant  flow  rates  at  the  injector  orifices  can  be 
obtained#  Design  of  the  unit  was  based  directly  on  experience  with  the 
monopropellant  flow  modulator,  and  incorporates  several  improvements  ad- 
ditional to  the  use  of  pressure-lubricated  roller  bearings  on  the  connect- 
ing rods#  A 5 HP  U.  So  Motors  Vari-drive  has  been  procured  to  drive  the 
unit  at  speeds  from  3,000  to  12,000  RfW# 

ihe  thrust  stand  and  thrust-measurement  devices  are  identical  to  those 
used  in  the  monopropellant  system#  All  piping,  fittings,  valving,  etc#  is 
of  typo  30U  or  316  stainless  steel  in  order  to  provide  for  general  usage 
with  nearly  all  of  the  common  rocket  propellants.  *‘hen  using  liquid  oxygen, 
or  any  other  low-temperature  propellant,  provision  has  been  made  for  pre- 
oooling  the  entire  feed  line  right  up  to  the  rocket  motor  injector  in  order 
to  avoid  the  hard  starts  somstimea  occasioned  by  vapor  in  the  feed  line# 

A>he  firing  diagram  fear  the  bipropellant  motor  is  shown  in  Figure  3o 
Special  features  include  pyrotechnic  igniters  utilizing  a safety  fuse  to 
prevent  propellant  valve  operation  in  case  of  igniter  misfiring,  a safety 
chamber  pressure  switch  which  will  shut  off  the  propellant  valves  la  case 
of  a rocket  motor  misfire,  and  automatic  pressurization  of  fuel  and  oxidant 
tanks  to  a preset  value.  Mo  special  provision  is  necessary  to  obtain  oxi- 
dant "la ad"  at  start  or  "lag"  at  shutdown  on  the  liquid-oxygen  system, 
sinoe  oxidant  is  automatically  introduced  into  the  rocket  motor  at  these 
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HI  APPARATUS  (Cant’d*) 

B.  Bipropellant  Rocket  Motor  (Cont’d#) 

times  via  the  precooling  arrangement  mentioned  earlier# 

Figure  U illustrates  the  initial  rocket  motor  design 0 Ihe  chamber  is 
of  1-^-inch  thiok  copper,  tliree  inches  in  inner  diameter  and  lour  inches 
longo  It  will  Operate  uncooled  for  the  first  aeries  of  tests,  in  order  to 
facilitate  convenient  instrumentation  placements  The  nozzle  is  also  of 
copper,  but  is  provided  with  water  passages  designed  for  cooling  at  uniform 
heat  transfer  rates#  The  injector  is  of  tne  one-on-one  impinging  jet  type, 
using  12  pairs  of  holes  located  on  a concentric  circle.,  -izes  of  tnc  fuel 
and  oxidant  orifices  are  identical,  the  proper  mixture  ratio  being  obtained 
by  variation  in  feed  pressures.  The  basic  design  of  this  injector  is  pat- 
terned after  a model  used  extensively  by  Reaction  motors,  inc#  for  perfonr*- 
ance  evaluation  and  instability  studies# 

Additional  instrumentation  components  consisting  of  a ^-channel 
electronic  Tube  Corporation  oscilloscope,  a U-channel  modification  of  the 
present  2-channel  Anpex  tape  recorder,  and  additional  pressure  pickups 
*nd  intermediate  equipment  have  beer,  orderod.  The  basil  instrumentation 
scheme  will  be  similar  to  that  used  on  tne  monopropeiiar.t  motor,  with  the 
exception  of  trie  heat  transfer  instrumentation  for  use  on  later  bipropel- 
lant test  operations# 

C#  Instrumentation 
1#  D.C#  Amplifiers 

The  problem  of  obtaining  D#  C.  amplification  conforming  to  project 
requirements  has  been  only  partially  solved#  T'nrec  channels  of  a chopper- 
feedback  design  have  been  developed  and  tested  and  ap.n-ar  to  ee  satiafact- 
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in  APFAttATUS  (Cont'do) 

C,  Instrumentation  (Cont'do) 

1*  Dq  C»  amplifiers  (Cont'do)  ' 

ory  from  linourity*  drifts  and  differential  amplification  standpoints!  but 
pnf  f cfp  from  an  excessive  noise  level*  Anese  instruments  were  used  to 
record  the  D*  C0  pressure  levels  described  later,  but  axe  apparently  un- 
suited for  use  with  low-level  A.  C,  signals  such  as  are  characteristic 
of  the  pressure  pcikup  -during  runs  with  modulated  flow  rate#  Tne  pres- 
sure-lag data  were  thus  recorded  without  benefit  of  D#  C,  differential 
simplification* 

It  is  believed  that  the  high  noise  level  is  not  an  inherent  feature 
of  this  design  of  amplifier!  but  rather  is  a function  of  the  detailed 
circuit  development*  During  the  last  month!  two  cownarTAal  amplifiers 
which  use  the  chopper-feedback  principle  and  which  aro  chained  to  be 
suitable  for  this  application  have  appeared  on  trio  market!  and  wilt  be 
investigated  as  possible  solutions  to  the  problem.  Meanwhile*  trie 
instrumentation  system  of  Figure  'S  was  calibrated  and  used,  to  record  the 
pressure  time- lag  data* 

2*  Overall  System  for  Noraopropeli?jnt  Treasure  Time  lag  measurement 

Figure  5 shows  the  instrumentation  usod  to  measure  the  time  lag 
between  a change  in  injector  pressure  and  the  resulting  variation  in 
chamber  pressure*  Two  differential  pressure  piclcups  were  used,  one  mount- 
ed in  one  of  the  injectors,  and  one  in  the  chamber  just  adjacent  to  that 

t 

injector*  The  signal  from  each  pickup  T.jas  then  split  four  ways*  One  branch 


• ■•*'*«***  , 


3*#  U*  Auro*  Lngg#-  Upt#  No#  226o  36 

111  APPARATUS  (Cont’d.) 

| C#  Instrumentation  (Cont’d#) 

2#  CVerall  Syat«n  for  rionopropelJ ant  treasure  lime  Lag  Measurement  (Cont'da) 

lod  directly  to  the  four-channel  oscilloscope,  upon  one  channels  of  which  a 
200  cpa  timing  trace  was  impressed  by  an  audio  oscillator  calibrated  with 
a Berkeley  precision  electronic  counter#  The  second  signal  branch  led  to 
a cathode  follower  and  thence  to  a Leeds  and  Northrup  recording  potent!- 

\ 

j ometer#  The  third  branch  wa3  passed  through  the  chopper-stabilised  D#  C# 
amplifier  and  wa3  recorded  on  the  Hathaway  oscillograph iJ  The  final  branch 
of  th  3 signal  from  each  pickup  was  amplified  by  a constant-phase  A.  C* 
amplifier  and  rooorded  on  magnetic  tape# 


The  input  impedance  of  each  of  the  instruments  described  abovo  was 
well  into  the  megohm  range,  and  hence  feedback  between  instruments  was 
negligible  compared  with  the  output  from  the  2!? 00-ohm  pickup#  The  L,*odo 
and  Northrup  potentiometer  records  were  used  to  obtain  the  D#  C#  level  of 
the  two  pressures,  viiiie  the  0#  C0  Hathaway  oscillograph  trace  was  used 
primarily  to  provide  a record  of  run  history  arid  to  show  transients  or 
other  unlooked-for  events#  Accurate  tine  values  were  taken  from  the  A#  C# 
Ampex  tape  recordings,  and  were  checked  with  oscilloscope  readings  (as 
reoorded  by  tho  General  Radio  atrip-film  camera)#  The  oscilloscope  and 
inpex  recordndldo  supplied  pressure  amplitude  measurements.  Keise 
Ibourdon-tube  gauge  readings  v-ere  used  to  check  the  steady-state  componj* 
ant  of  each  pressure  recording# 

i 

~ “ 3s  Pressure  Pickups 

•3 

Three  differential  pressure  pickups  incorporating  the  Improved  coolant 


j III  APPARATUS  (Cont'd.) 

0*  Instrumentation  (Cont'd*) 

i 

3s  Pressure  Pickups  (Cont'do) 

i 

] 

water  passage  design*  extended  Iin*ar  range*  and  back-pressure  chamber 
sealing  modiiications  have  been  delivered*  and  have  been  used  to  obtain  all 

j 

tost  data  included  in  this  report*  Calibration  tests  show  what  first  ap- 
I peered  to  be  a large  amount  of  hysteresis*  but  which  was  demonstrated  to 

I oocur  only  when  high  back  pressures  were  applied  without  any  test  pressure* 

] 

Calibrations  mads  under  simulated  operating  conditions  were  entirely  satis- 

i 

factory*  as  will  be  illustrated  later  in  the  report*  Pour  additional 
pickups  of  this  type  have  been  ordered*  and  a slight  modification  in  dia- 

t 

phragn  spacer  design  is  expected  to  eliminate  even  the  apparent  hysteresis 
Wbioh  occurs  under  the  artificial  calibration  condition  just  described* 

Thermal  drift  and  "screaming"  heat  transfer  motor  tests  on  an  absolute 
pressure  pickup  with  the  improved  coolant' passage  design  were  made  at  our 
request  by  the  NACA  Lewis  Alight  Propulsion  Laboratory*  A differential 
piokup  was  also  sent  to  the  Bell  Aircraft  Corporation  for  evaluation  of 

: thermal  drift*  vibration  sensitivity*  and  rug^edness  under  typical  rocket 

i 

actor  test  conditions.  Aie  NACA  data  are  included  later  in  the  present 
report*  while  the  Bell  tests  will  be  reported  upon  when  complete* 

it*  Flownetere 

i 

%e  Mittelmann  electromagnetic  flometer  was  operated  for  an  extended 
' period  without  achiering  so  much  as  a steady-state  calibration*  Multitudes 
of  olrouitry  difficulties  ranging  from  poor  Contacts  to  resonance  in  the 
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XU  APPARATUS  (Cont'd.) 

C,  Instrumentation  (Cont'd* ) 

i 

U«  Flowmeters  (Cont'd.) 

Cables  wore  encountered,  and  the  instrument  is  to  be  returned  to  the 
vendor  far  revision  and  repair.  Due  to  its  inability  to  measure  even  a 
steady  flow  rate,  no  evaluation  of  the  flowmeter's  dynamic  performance 
could  be  made,  and  the  only  recommendation  that  can  be  forwarded  at  this 
tine  is  that  a good  deal  of  electronic  circuit  cleanup  is  necessary  be- 
fore the  meter  beoomes  suitable  for  use  in  any  flow  measurement  applica- 
tion* 

The  il  mass  f lowneter  has  been  received  after  revision  by  the  vendor 
to  correct  its  excessive  leakage*  Preliminary  checks  indicate  its  suit- 
ability for  use  in  calibration  tests,  but  the  leakage  rate,  although  much 
smaller,  still  precludes  the  instrument's  use  on  direct  rocket  motor 
tests*  *’ull  dynamic  water  calibrations  in  conjunction  with  the  hot-wire 
liquid  flow  phase-meter  discussed  in  the  previous  wuartarJy  Report  will 
be  made  on  the  monopropellant  flow  stand  at  the  conclusion  of  t*he  present 
series  of  rocket  motor  tests* 


* « J • . \ tV  O o 
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iv.  r:i‘’0'ti"noiv  a;.d  data 

A«  j ionopropellant  Rocket  Tests 


'ftie  data  taken  during  the  current  aeries  of  toot  runs  are  not  yet 
considered  to  be  of  sufficient  accuracy  for  comparison  with  theory,  Al- 
though  no  direct  measure  of  instantaneous  flow  is  being  taken  on  these 
test35  it  is  ejected  that  subsequent  theoretical  corrections  checked  by 
flow-stand  calibrations  vri.ll  permit  use  of  the  instantaneous  pressure 
drop  readings  for  flowrate  d^erminations , following  the  theory  outlined 
in  Appendix  A of  the  Third  Quarterly  Progress  Report# 

Tho  difficulty  of  obtaining  a direct  rcjcord  of  the  instantaneous 
« 

injector  pressure  drop  through  a suitable  D#  C#  amplifier  makes  it  prefer- 
able to  obtain  the  corresponding  data  analytically  from  more  elementary 

\ 

measurements#  In  fact,  determination  of  the  injector  pressure  drop  can 
be  made  using  only  the  moan  levels,  amplitudes  of  oscillation,  and  phase 
difference  between  injector  pressure  and  chamber  pressure#  If  the  oscil- 
lating portions  of  the  injector  and  chamber  pressures  are  sinusoidal  with 
respect  to  time,  the  injector  pressure  P^  may  be  written 

_ *s/ 

=.  Pc  PL  5) n 03 

The  chamber  pressure  P then  becomes 

c 

5> * v-\  (U3i  1 oc) 

where  o(.  58  measured  phase  difference  between  injector  and  chamber 
pressure  oscillations 

P P.  = mean  values  of  injector  and  chamber  pressures  respectively 
^ sa  amplitudes  of  oscillation  of  Injector  and  chamber  pressures. 


respectively# 
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IV.  INFORMATION  AND  DATA  (Cont‘d8) 

A,  Monopropellant  Rocket  Tests  ( Cont ' d, ) 

The  instantaneous  pressure  drop  AP  across  the  injector  may  be  written 
in  the  form 

= pL_  pc  = AP  1-  £?  t *) 

where  = phase  difference  between  ?4  and  Ap. 

But,  we  are  ultimately  interested  in  determining  the  pressure  perturbation 
as  a function  of  the  mass  flow  perturbation.  Assuming  for  the  moment  that 
the  relationship  butween  flow  rate  and  AP  were  known,  this  means  that  we 
meed  the  relation  between  P0  'and  the  quantity 


Av>  - AP 

& 


Hence,  we  are  intere  sted  in  the  ratio 


This  can  be  written  in  the  form 


where  the  amplitude  of  R of  the  required  ratio  is  given  by 
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IV.  IHDTBUM^XATI'Jl'i  AND  DATA  (Cont»d.) 

A«  Honoprop allan t Kockot  Te3t3  (Cont  >d0 ) 


and  the  phase  difference  ^ betweon  Pc  and  AP  is 


& 


Co&  <*. 


- cMO 


'Thus,  the  necessary  relationship  between  chamocr  pressure 
perturbation  and  injector  pressure  drop  perturbation  (wuich  can  oe 
correlated  to  the  flow  perturbation)  is  obtained  directly  from  tne 
mean  values,  amplitudes  of  oscillations,  and  phase  difference  between 
injector  and  chamber  pressure.  This  determination  can  bo  made  with- 
out the  use  of  a hijh-response  D.  C.  amplifier,  requiring  onlv  that 
the  pressure  oscillations  be  sinusoidal  in  nature. 

Typical  oscilloscope  recordings  of  injector  ana  enamour  pressur"  -jp, 
are  shown  in  figure  6.  The  uppermost  trace  is  a 200  cps  timing  trace, 
the  middle  one  is  injector  pressure,  and  the  lower  one  is  cnamoer 
pressure. 

forthcoming  cold  flow  stand  tests  wilu.  oe  used  to  determine  whether 
the  hi-jh-freqnency  hash  (about  1,000  cps ) which  appears  on  tno  injector 
pressure  trace  is  due  to  coiobustion  or  to  fluid  line  osciiuations.  In 
any  case,  before  entirely  adequate  phase  and  amplitude  measurements  can 
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IV  INFURMATIJL1  AMD  DATA  (Cont‘d.) 

A0  honopropeliant  Socket  Tea  Ay  (Cont'd,) 

be  drawn  from  this  data,  it  will  be  ne<?essi\ry  to  remove  the  hash  by  means 
of  a suitable  low-pass  Tiltar0 

The  excessively  high  amplitude  of  ti.e  injector  pressure  oscillation  in 
the  current  tests  results  from  the  fact  that  the  largoct-sizo  piston  is  being 
used  in  the  f lxw  modulating  unit,  iaoB,  the  one  designed  for  900  psi  chamber 
pressure  operation,  Ahis  one  piston  has  beon  used  at  all  chamber  pressures 
(the  run  illustrated  in  Figure  6 is  at  300  psi)  only  as  a natter  of  - conven- 
ience, but  "hen  the  final  sories  of  tests  is  made,  the  proper  piston  size 
for  each  value  of  the  .mean  propellant  flowrate  will  be  u3ed. 

13.  Freagure  Pickups 

Results  of  teals  made  at  tmo  '.AC A on  heat- transfer  and  temperature 
characteristics  of  the  mi- Liu  double-diaphragm  pickups  are  included  in 
Figures  7,  0,  and  9o  Figures  7 and  B snow  tine  coolant  water  flow,  coolant 
temperature,  and  heat  transfer  rate  for  a non- screaming  and  screaming  run, 
respectively,  Ahe  data  of  Figure  B are  incomplete  due  to  burnout  of  tin* 
rocket  motor  during  screaking.  lio  damage  to  tea  viekup  war.  observed.  ( Ihe 
heat  transf  r rate  is  based  on  the  overall  exposed  area  of  t e pickup.  If 
only  the  diaphragm  area  is  considered,  tnc  heat  transfer  rates  shown  in  the 
figures  should  be  multiplied  by  about  2,2,  The  actual  heat  transfer  rate 
probably  falls  somewhere  between  the  tvo. ) 

Figure  9 shows  the  tespearoture  performance  of  an  absolute  pickup  up  to 


IV.  L.FiAL- iATIuli  AMD  D.a'A  (Coat'd.) 

B.  Pressure  Pickups  (Cont  'd«  ) 

500CF.  Aiese  data  are  gone  /hat  questionable,  due  to  outh  the  long-term 
exposure  to  furnace  temperature  and  the  low-accuracy  pressure  gauges  used 
on  the  tests*  but  illustrate  the  general  nature  of  the  pickup's  performance. 
Actually,  it  is  riot  possible  to  detr>rnune  from  these  data  whether  the  drift 
is  due  to  temperature  or  change  in  sensitivity,  and  additional  temperature 
tests  will  be  made  both  to  clarify  this  point  and  to  supply  a more  satis-  * 
factory  temperature  characteristic,  during  the  temperature  drift  tests, 
the  signal  dropped  suddenly  to  zaro,  and  the  pickup  (absolute  pickup  ;/h ) was 
returned  to  the  vendor  for  internal  examination.  Mo  exterior  damage  was  ob- 
served. Another  live  pickup  (#9)  was  3ent  to  NACA  for  evaluation  on  a 
screaming  rocket  motor. 

AU  available  ini' or;  vat  ion  on  this  pickup  design  vtilu.  be  collected 
and  brought  up  to  date  for  presentation  at  the  becond  Symposium  on  liquid 
Propellant  Rocket  Combustion  Instability  and  will  appear  as  an  appendix  to 
the  next  Quarterly  Report 0 

Co  Theory 

A theoretical  analysis  of  the  effect  of  nozzle  impedance  on  low-frequency 
combustion  instability  in  liquid-propellant  rockets  is  included  as  Appendix  A 
to  the  present  report#  This  analysis  indicates  that  the  low-frequency  stabil- 
ity of  a rocket  motor  is  improved  by  decreasing  the  nozzle  angle  of  converg- 
ence, irrespective  of  wuother  this  is  done  by  decreasing  the  contraction  ratio 


O 
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IV.  INPcRI-iATlOfi  AND  DiiTA  (Cont'd.) 

C*  Theory  (Cont'd0) 

or  increasing  the  length  of  thoconvergant  section. 

V.  DISCUSSION 

Detailed  discussion  of  the  monopropellant  test  data  will  be  postponed 
until  analysis  of  the  current  series  of  runs  has  been  conpleted.  Evaluation 
of  results  to  date  on  the  pressure  pickup  development  program  will  appear  in 
the  next  quarterly* 
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SCHEMATIC:  BIPROPELLANT  FEED  SYSTEM 


TEFLON  CUP  PACKING 


FIGURE  2: 

BIPROPELLANT  FLOW  MODULATING  UNIT 


PRESSURE  PICKUP  Hm  f 


bipropellant  ROCKET  MOTOR  N2  I 
(LIQUID  OXYGEN-ALCOHOL) 


»7|UJ 
W(L  ^ 

Soq: 

goo 

xgt 


0 — 5 

>ot 

rco 

»~o 


M 

P~ 

1— 

, CL 

UJ  « 
p»  Q 

5 or 

Is 

Su 
^ o- 


tf) 

i o: 

H W 

Zuj  t 

5cn^ 

f~  < Q_ 

5i5 
Z cl  ^ 
O ^ 
O o 
C 


CO 

cr 

LU 
Oh- 
Zuj 
Q5 
cr  o 
op 
05 

uj.f; 
a:  w 


co 

UJ* 

QUJ 

O £ 
X o 

h-  -I 
< _J 

O o 


a 

ocx 

O Q_ 

o< 

□:  ^r) 

.O 


z o 
< co 
x o 

0 

1 

X 

CO 


CO 

ujct 

QUJ 

xo 

H-J 

<_j 

o o 

u_ 


Q?r 


cr  uj 
O o 
h-  Z 
< UJ 
_J  X 


or 

UJLlI 

CD 

Zg- 

Z 

zO 

<o 

x<3 

h- 

uO 

1 -1 
cr  _J 

26 
0 CO 

U-O 

U_  UJ 
' 5 

Q=« 

cr  o 


cr£ 

UJ* 

£w 

geo 

X ct; 
<->£ 


Uiuj 

OrZQ 

=>ujco 

XCCCO 

OUJUJ 

K^-CC 

UJCL 

o: 


/LI-LIU  \ POWER  SUPPLY 

/ PRESSURE 
/ PICKUP 

FIGURE  5 

SCHEMATIC ■ MONOPROPELLANT  PRESSURE -LAG  INSTRUMENTATION 


COOLING  WATER  AT  (°F) 


HEAT  TRANSFER  TO  LI -LIU  DUMMY 
PRESSURE  PICKUP  (NON- SCREAMING) 
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HEAT  TRANSFER  RATE , Q (BTU/m  - sec.) 


COOLING  WATER  AT  (°F) 


HEAT  TRANSFER  TO  LI -LIU  DUMMY 
PRESSURE  PICKUP  (SCREAMING) 
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FIGURE  8 TIME  (SECS) 


THERMAL  DRIFT  OF  LI-LIU  PICKUP  NO  9 
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PSNDIX  A 


"LCKV  FRKiUENCY  CGMBUSTIUN  STABILITY 
OF  U14JID  ROCKET  MOTOR  WITH  DIFFLRLiN'T 

NC2ZLLS" 
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1. 

LOW  FREQUENCY  GOMBUSTIGN  STABILITY  CP  LIQUID  ROCKET 
MOTOR  WITH  DIFFERENT  NOZZIiS 

I*  INTRODUCTION 

,ft*evious  analyses  of  low  frequency  oscillations  (1)  to  (?)*  are  made  on  the. 
assumption  of  quasi-steady  gas  flow  in  deLaval  nozzle  through  which  combustion 
gases  are  accelerated  and  ejected  at  supersonic  speeds*  The  boundary  condition 
as  deduced  from  this  quasi-steady  flow  assumption  is  that  the  Mach  number  of  the 
gas  flow  entering  the  de Laval  nozzle  is  constant  so  that  the  rate  of  flow  of 
burned  gas  out  of  the  combustion  chamber  into  the  deLaval  nozzle  is  directly 
proportional  to  the  stagnation  pressure  and  inversely  proportional  to  the  square 
root  of  the  stagnation  temperature  of  the  ourned  gas  entering  the  nozzle*  Taien 
(0)  questioned  the  validity  of  this  assumption  and  analyzed  the  transfer  function 
of  the  nozzle  for  several  special  cases*  Crocco  (9)  extended  Tsien's  treatment 
with  linear  steady  state  velocity  distribution  in  subsonic  portion  of  the  nozzle 
and  determined  the  specific  acoustical  admittance  ratio  of  the  nozzle  flow  for 
considerable  range  of  the  reduced  frequency  parameter  p • The  results  show  that 
for  small  values  of  p the  real  part  of  the  admittance  ratio  is  essentially  constant 
and  is  equal  to  the  vaiue  corresponding  to  that  of  quasi-steady  nozzle  flow* 
However*  the  imaginary  part  of  the  admittance  ratio  is  directly  proportional  to 
the  reduced  frequency  p , and  the  rate  of  increase  of  the  imaginary  part  with  p is 
quite  considerable  especially  when  the  entering  Mach  number  is  small  as  is  for  the 
practical  case*  Furthermore,  the  magnitude  of  p for  chugging  frequencies  with 
a conventional  nozzle  is  not  really  very  small  but  is  of  the  order  of  one  tenth 
and  can  be  considerably  bigger  if  the  deLaval  nozzle  is  exceptionally  long.  The 
object  of  the  present  paper  is  to  investigate  the  effect  of  the  deviation  of  the 
nozzle  flow  from  being  quasi-steady  on  the  stability  of  ldw  frequency  oscillation 
on  liquid  propellant  rockets 0 


* Humber  in  parenthesis  refers  to  references* 
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II.  Fundamental  Equations 


The  status  in  the  combustion  chamber  of  a liquid  propellant  rocket  is 
extremely  complicated.  For  simplicity,  we  shall  first  consider  the  monopropellant 
case  and  adopt  the  following  assumptions  in  accordance  with  previous  authors  without 
further  discussion. 

a«.  The  gas  pressure  inside  the  combustion  chamber  is  practically  uniform 
everywhere  in  the  combustion  chamber  and  at  any  instant.  In  unsteady  state 
operation,  the  chamber  pressure  oscillates  about  the  mean  or  steady  state  value 
as  a whole. 

b.  The  adiabatic  flame  temperature  of  the  burned  gas  is  assumed  to  be 
constant  and  independent  of  the  small  variations  of  pressure  under  which  combustion 
takas  place.  The  temperature  of  the  burned  gas  is  further  assumed  to  remain 
unchanged  throughout  its  motion  in  the  combustion  chamber  irrespective  of  the 
chamber  pressure  oscillation. 

c.  The  time  lag  of  all  the  propellant  elements  are  assumed  to  be  the  same. 

The  first  two  assumptions  enable  us  to  neglect  the  equations  of  conservation  of 

momentum  and  of  energy  for  the  flow  of  the  burned  gas  in  the  combustion  chamber. 

The  dynamics  of  the  burned  gas  flow  is  governed  only  by  the  equation  of  conservation 

of  mass,  that  is,  the  rate  of  burned  gas  generation  Ap  (t)  must  be  equal  to  the 

sum  of  the  rate  of  ejection  Ae  (t)  of  burned  gas  out  of  the  combustion  chamber  into 

cl  I*\%. 

the  deLaval  nozzle  and  the  rate  of  accumulation  5c"  of  burned  gas  in  the 
combustion  chamber.  All  notations  are  adapted  from  Crocco  (k)  except  where  new 
symbols  are  necessary.  The  equation  of  mass  balance  is  reduced  to  the  following 
dimensionless  forms 


J1  f/t£  (g)  ,/u.i  (g) 


t_ 

where  ^ ? s ^ is  the  reduced  time  with  <9^  3 y 

the  average  gas  residence  time  in  the  combustion  chamber. 


-Jl±L 

>-/ 


L 

C 


(2.1) 


3. 


We  shall  assume  that  the  entire  time  lag  is  pressure  sensitive  and  the 
dependence  of  "i V on  pressure  is  defined  as 

? j=>n'  (t‘)  d.t  f - 

t-T 

where  n is  the  index  of  interaction  in  terms  of  instantaneous  pressure  including 
the  effect  of  temperature  variations  of  the  burned  gas*  It  has  been  shown  ( k ) 
that 

(2.2) 

The  fractional  variation  of  the  mass  flow  rate  out  of  the  combustion  chamber 
depends  on  the  transfer  function  Nn  of  the  deLaval  nozzle  with  supersonic  exit 
velocity*  For  quasi-steady  nozzle  flow  and  isentropic  oscillations  of  the  burned 

gas*  we  have  ( 

_ A'  ~7~  ‘ _ y + ? op 

j^C.  - f°  ^ ~r  ~ A y 

with  the  nozzle  transfer  function  based  upon  fractional  pressure  variation  given 
as  £-^2-  which  is  equal  to  unity  when  the  temperature  variation  of  the  gas  is 

z-  y 

neglected*  The  specific  admittance  ratio  is  defined  as  the  ratio  of  the 
fractional  variation  of  gas  velocity  to  that  of  gas  density  at  the  entrance  of 
the  defeval  nozzle  and  can  be  represented  for  low  frequency  oscillations  as 


U.K 


(2.3) 


where  is  the  dimensional  angular  frequency  of  neutral  oscillation  and  ux  is 
the  constant  velocity  gradient  c*  the  gas  flow  in  the  subsonic  part  of  the  nozzle* 
The  constant  k for  isentropic  oscillations  can  be  obtained  from  reference  (9) 


u» 


as 

(*.u 

in  which  •¥  can  be  neglected  in  accordance  with  assumption  b/.  The  quantity 
£ & is  defined  as 

with  Alt  representing  the  Mach  number  of  the  burned  gas  entering  the  nozzle*  The 
quantity  ? ^ is  usually  less  than  0*1  for  conventional  liquid  rockets*  The 
magnitude  of  for  a deLaval  nozzle  with  converging  portion  of  1 foot  long 

and  for  an  oscillation  of  100  cycles  per  second  is  of  the  order  of  0*2  for 
example.  Thus  we  see  that  the  imaginary  part  of  the  specific  admittance  ratio©(^ 
is  quite  significant  as  compared  to  its  real  part*  The  transfer  function  of  the 
mass  flow  based  on  fractional  pressure  variation  at  the  entrance  is 

K - i ' * * * u* 

which  reduces  to 

* i + M-*-  4“  = b 00 

idien  Y = 1 under  the  assumption  of  negligible  temperature  variation.  The 
proportionality  constant  b is  dimensionless  and  is  equal  to  /Cog.  ■ & x)  • 

When  the  steady  state  velocity  in  the  subsonic  portion  of  the  nozzle  is  linear f 
the  value  of  k is  given  by  equation  (2«ii)  and  the  factor  b is  given  as 

) - _ _ //  Jb*'  Ejm.  1 _ / 

/-f£  • aK  ' (2.5) 

It  is  conjectured  that  for  practical  nozzle  shape  where  the  velocity  distribution 


'+-!£L 


is  not  linear,  this  constant  b will  still  be  independent  of  the  frequency  of 
oscillation  and  of  the  same  order  of  magnitude  as  that  with  linear  velocity 
distribution  in  nozzle#  It  is  thus  seen  that  the  absolute  magnitude  of  the 
transfer  function  will  ndt  be  significantly  different  from  unity,  which  is  the 
quasi-steady  value,  but  the  phase  lead  of  the  mass  flow  oscillation  is  quite 
considerable  even  for  low  frequency  oscillations* 


ol 


If  we  write  the  nozzle  transfer  function  in  operational  form  - t + b 
while  investigating  the  stability  of  the  solutions  of  exponential  type,  then 

--  Cl  + b 3*8  ) <-lfh  ft  ) t 


A 


(2.6) 


where  ^ is  replaced  by  9 since  we  assumed  that  the  fractional  variation  of 
pressure  at  any  instant  is  the  same  everywhere  in  the  combustion  chamber* 
Introducing  equations  (2,2)  and  (2,6)  into  equation  (2.1),  we  obtained  the 
equation  of  mass  balance  as* 

- 0 tii*-)  f-f* 

Accumulation  Ejection  Generation 


(2.7) 


The  formulation  for  the  analysis  of  the  oscillation  problem  will  require  a 
relation  between  the  fractional  variation  of  the  injection  rate  and  the 
fractional  variation  of  chamber  pressure*  This  relation  depends  upon  the 
dynamic  characteristics  of  the  feeding  system  in  response  to  the  chamber  pressure 
variation*  The  dynamics  of  the  feeding  system  will  be  formulated  following 
Crocco  (U)  and  Tsien  ($),  with  the  following  assumptions* 

(a) .  The  pressure  drop  due  to  wall  friction  in  feed  line  is  negligibly  small 
compared  to  the  overall  drop^of  the  feeding  system. 

(b) .  The  effect  of  the  elasticity  of  the  feeding  line  can  be  represented  by 
an  equivalent  concentrated  spring  capacitance  <2$  located  at  a distance  4 . JL 
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downstream  of  the  pump  delivery,  where  JL  is  the  equivalent  lengtn  of  the  feeding 
line.  The  spring  constant  of  the  capacitance  is  defined  as  the  change  in 
volume  of  the  feed  line  produced  by  unit  pressure  rise  over  the  entire  feed  line. 

(c) .  The  feed  pump  responds  to  oscillations  of  delivery  pressure  in- 
stantaneously and  the  fractional  variation  of  mass  delivery  rate  is  directly 
proportional  to  the  fractional  variation  of  the  delivery  pressure  with  the 
proportionality  constant  D defined  as 

— ■ = - d -4~- 

P * (2.8) 

and  D is  assumed  real  and  therefore  can  be  evauated  from  the  steady  state 
performance  curve  of  the  pump.  v 

(d) ,  A control  capacitance  Cc  excited  by  a feed  back  circuit  is  introduced 
just  upstream  of  the  injector.  The  characteristics  of  the  feed  back  circuit  is 
described  by  the  transfer  function  F ( ~oLW~ / ** 

f (a*-)  7 " srSi 


(2.9) 


A schematic  diagram  of  the  liquid  rocket  is  given  in  figure  1.  A straight  forward 
process  of  elimination  leads  to  the  following  dimensionless  equation  of  the 
feeding  system  dynamics  relating  and  as 

[i+pe  j Efr  ]+[d(p+&) Zs ^ (f-p)cr+& ) 

s t* Ey.  O-p  )zt?]  f (ZjT  )J  J 

. ,* 

-h  £ li+x>  Cp*S )] *■[&£  (-Pf'/x)s zj  + \j>i e ( /-yXr+Xi) (2ol0) 

+ = o. 

p - 

yhere  ? » JAp  is  the  pressure  drop  parameter,  a relative  measure  of  the 

ad  f=a  yC 

pressure  drop  across  the  injector.  S » is  the  elasticity 


TAK1U,  LINE 

TANK  CAPACITANCE  SERVO 


NOZZLE 
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parameter,  a ratio  of  the  rate  of  mass  accumulation  in  the  line  capacitance  due 

to  a rate  of  pressure  increase  of  / £>£.  to  the  mean  mass  flow  rate  of  the 

JL  ■>»i 

system,  J = 3 a jb  ■ At-  is  the  inertia  parameter,  a ratio  of  the  time  required 
to  accelerate  a given  element  from  rest  to  the  state  of  motion  in  the  feed  line 
under  the  pressure  drop  of  3 A jb  as  compared  to  the  gas  residence  time 

Equation  (£.10)  essentially  defines  the  ratio  f*- j.  / ? in  operational  form 
which  is  the  transfer  function  of  the  entire  feeding  system  iii  response  to  the 
chamber  pressure  variation.  Substituting  from  equation  (2.10)  into  (2.7), 

ws  obtain  a single  hystero-ordinary  differential  equation  with  constant 
coefficients  for  the  fractional  pressure  perturbation  f as  a function  of  the 
reduced  time  H « The  characteristic  equation  for  the  investigation  of  the 
stability  of  the  solutions  of  exponential  type  is  obtained  by  replacing  the 
operator  sdh  with  the  characteristic  value  s.  It  is  obvious  that  the 
dynamics  of  the  feeding  system  which  creates  the  variation  of  the  injection 
rate  is  responsible  for  most  of  the  algebraic  complications  for  the  solution  of 
the  problem* 
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III.  TRANSFORMATION  OF  CHARACTERISTIC  EQUATIONS  WITH  DIFFERENT 

ROCKET  NOZZLES 

% 

It  is  interesting  to  note  that  the  correction  factor  b for  the  phase  lead 
component  of  the  nozzle  transfer  function  which  represents  the  deviation  of  the 
nozzle  flow  from  being  quasi-steady  can  be  absorbed  by  defining  a new  characteristic 
time  instead  of  the  average  gas  residence  time  ^ as  defined  by 

Crocco#  The  dimensionless  parameters  which  involve  the  characteristic  time  should 

of  course  be  correspondingly  modified#  Thus 

t/^/  * ?/  0+&) 

£ - - S / C'+’b) 

['  --  a#  = c/a  + b) 

~y~  - yt  -rrt  / £L  A ■ A ' <5^  - 3" / C I + & ) 

t ('A7')  • &y  = r (di  ')/ci+6 ) 

and  ou'  = - (_/ tb  ) (3.1) 

Those  dimensionless  parameters  P,  D,  and  n which  do  not  involve  characteristic 
time  remain  unchanged. 

With  this  similarity  transformation,  equations  (2.7)  and  (2.10)  can  be  easily 
verified  to  be  completely  independent  of  b and  are  identical  with  the  form  of 
these  equations  with  b = o corresponding  to  a "very  short  nozzle"  where  the 
quasi-steady  nozzle  flow  assumption  is  approximately  valid.  For  different  delaval 
nozzles  the  values  of  b will  be  different,  but  its  effect  can  be  accounted  for 
simply  by  changing  the  characteristic  time  according  to  ® (1+b).  The  solution 

of  the  characteristic  equation  remains  identical  in  terms  of  the  dimensionless 
quantities  defined  in  equations  (3.1).  Therefore  the  solutions  and  the  stability 
boundaries  given  by  previous  authors  for  the  case  b = 0 in  terms  of  these 
dimensionless  quantities  represent  the  similar  or  the  universal  solutions  for 
rockets  with  given  feeding  system  and  combustion  chamber  but  different  delaval 
nozzles. 


?. 

The  characteristic  equation  becomes, 

(s7/-*.*ne”*/*;  £[t*d  (r+& ')]^[DE  (?&)+’*]*' 

t [dtW  Ci-^CP+'/a.)*  j'e  '%]  s'Vx'V^  0-£)s/3J 

= P[/+DE'(r+  s,aJ 


+[*>  (P+  £)+j'i  v 2)3/£/  0-2.)CP+  %L)s'a+T'*E/2.  (/-#)*'*] s'f'  (s)  • (3>  2 J 

Prom  equation  (3*2),  we  observe  that  the  dimensionless  quantities  E/and  j' occur 
only  as  E/  and  J x <su^  when  s'  - ji>  for  the  determination  of  the 
stability  boundary*  Equations  (3*1)  show  that  these  products  E ' M ' and  J ' ^->'are 
independent  of  b and  equal  to  E ^ and  J cu  respectively.  The  quantity  <*>'  S 
which  represents  2 7T  times  the  ratio  of  time  lag  to  the  period  of  oscillation  is 
also  independent  b with  S * - <4*.  £.*  • Consequently,  the  effect  of  non- 
zero values  of  b enters  only  through  s'  and  sVYs'Jsi,  f'  ) 

in  equation  (3.2).  The  critical  value  of  ^ for  neutral  oscillations  of 

frequency  is  a function  n,  P,  D,  y,  E / oo'#.  and  the  coefficients  of 

f'Cs')*  If  the  feeding  system  of  the  rocket  is  not  servo-controlled  by  feed  back 
link,  then  F/=  F = o*  The  critical  values  of  and  co^_  are  therefore 

functions  of  dimensionless  quantities  which  are  all  independent  of  b.  for  such 
uncontrolled  systems,  the  effect  of  the  phase  lead  component  JLb^  can  be  directly 
obtained  without  determining  the  similar  or  the  universal  solution  of  the 
characteristic  equation.  Since  and  S*  are  constants  for  any  values  of 
b for  a given  uncontrolled  system,  and  since 
°°X  » G + b ) 


/ / 


>*  * a Tr 


Jl  • T-  - ^ ' 


C I + b ) 


that  for  increasing  b or  larger  phase  lead,  the  dimensional  frequency  </2. 


we  see 
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in  cycle a per  second  of  neutral  oscillation  is  decreased  and  the  dimensional 
critical  time  lag  is  increased.  It  can  therefore  be  concluded  that  for  rockets 
vith  feed  system  not  servo-controlled  the  phase  lead  component  of  the  nozzle 
transfer  function  tends  to  stabilize  the  system  towards  low  frequency  oscillations, 
and  larger  phase  lead  component  results  in  larger  stabilizing  effect  in  increasing 
the  critical  time  lag. 

For  rockets  with  feed  system  controlled  by  a given  feed  bac’:  3ervo-mechanism, 
the  coefficients  in  the  transfer  function  F(s)  are  different  from  those  in  F'Cs*). 
The  qualitative  effect  of  the  nozzle  constant  b cannot  be  drawn  without  further 
investigation. 
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IV.  BIPROPELLANT  ROCKETS 


In  bipropellant  rocket  systems  the  dynamic  behaviors  of  the  feeding  systems 
of  the  oxidizer  and  the  fuel  are  in  general  different.  The  variations  of  the 
injection  rates  of  oxidizer  and  of  fuels  in  response  to  the  same  chamber  pressure 
variation  are  different  with  a consequent  variation  in  the  oxidizer-fuel  or 
mixture  ratio  /-o  - ^ . Hie  adiabatic  flame  temperature  of  a given 

propellant  combination  depends  to  a certain  extent  on  the  mixture  ratio.  He 
stagnation  temperature  of  the  burned  gas  at  a given  position  in  the  combustion 
chamber  therefore  varies  with  time  as  the  chamber  pressure  oscillates  and  the 
stagnation  temperature  of  the  burned  gas  at  different  positions  in  the  chamber 
at  a given  instant  are  also  different.  His  variation  of  stagnation  temperature 
of  burned  gas , which  is  absent  in  monopropellant  rockets,  must  be  taken  into 
account  in  formulating  the  equation  of  mass  balance  in  the  combustion  chamber 
for  bipropellant  rockets. 


'A  ’Ai, 


0>.1) 


By  assuming  that  all  combustion  takes  place  practically  at  the  injector  end, 
and  that  all  propellant  elements  have  the  same  gas  residence  time  during  which 
all  these  particles  travel  from  the  injector  end  to  the  combustion  chamber  exit 
and  preserve  their  respective  temperature  at  the  instant  when  they  were  generated, 

. . ra+e 

Crocco  (a)  found  the  following  expressions  for  the  fractional  burning  perturbation 


*■  (*-S)  + n [ ?€*)-?  <>-£)] 


(4.2) 


with 


A*  * C£  + * ) f+0  + 


(4.3) 
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_/  n.~  / 

•ad  /' 

For  the  rate  of  mass  accumulation, 


2K  £a>  (z-S)-fr 


(b.k) 


For  the  fractional  variation  of  ejection  rate  evaluated  under  the  assumption  of 
° quasi-steady  nozzle  flow, 


vhere  2 K 


Ae  - f (^) - K [/S,  (z-s-/)-/'V  CHS-')] 


(^•5) 


represents  the  dimensionless  slope  of  the  adiabatic 


T drt. 

temperature  curve  for  different  mixture  ratio. 

Just  like  the  monopropellant  case,  the  fact  that  the  flow  of  burned  gas  in 
the  nozzle  is  not  quasi-steady  only  modifies  f*-  e.  , not  the  other  terms.  But 
owing  to  the  variation  of  the  stagnation  temperature  of  the  burned  gas,  there 
is  also  an  entropy  oscillation  of  the  burned  gas  entering  the  nozzle.  The  specific 
admittance  ratio  is  given  by  Crocco  (9)  as 


JkJt  - 


td- 

x o-*  u-x 


{h.6) 


for  a de Laval  nozzle  with  linear  steady  state  velocity  in  the  converging  portion 
vhere  0 and  O'  are  the  amplitudes  of  the  temperature  and  the  entropy  oscillations 
at  the  entrance  of  the  nozzle  respectively.  In  accordance  with  the  assumption  that 
the  temperature  oscillation  of  the  burned  gas  in  the  combustion  chamber  is 
neglected  when  the  oscillation  is  lsentropic,  # -1  is  neglected  as  compared  to 
unity.  Thus  the  first  approximation  of  phase  lead  component  of  the  admittance 

ratio  becomes  - 4 which  is  Identical  with  the  value  given 

/-  * c. 
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in  equation  (2.4)  for  monopropellant  case.  That  the  phase  lead  coniponent  is 
independent  of  the  ratio  wB  of  the  angplitude  of  entropy  to  temperature 
oscillation  is  correct  only  when  y ■ 1.  Since  V-/  for  practical  combustion 
gases  are  usually  0.2  or  0.3;  the  ratio  of  may  be  of  importance,  and  the 

magnitude  of  this  phase  lead  component  will  depend  on  many  thermodynamic 
properties  of  the  combustion  gases  even  for  an  Idealized  nozzle  with  linear 
steady  state  velocity  distribution. 

Hie  transfer  function  for  the  mass  flow  through  the  nozzle  is,  for  the 
bipropellant  case, 

~?f  ’ K * 0+«J  ( '-.-ft) 

Thus  in  accordance  with  the  assumption  of  / , we  have 


= C / ^ ^ d X ) C f e.  ) 


(4.8) 


/ • 


where  b «■  -1  -^=-  £Ttl,  vhich  is  the  same  for  monopropellant  case.  The 

* * c j-  K 

egression  for  --  % - ^ , when  b * 0 is  given  by  equation  (4.5)-  11“  the 

steady  state  gas  velocity  in  the  nozzle  is  not  linear,  we  would  of  course  expect 
a slight  change  in  the  expression  and  the  value  of  b.  The  equation  of  mass  balance 
becomes 

C/+6)  + 0-**)  f Cx-s) 

yu0  Cx-sy-C-zX-J:*"  Cx-$) 

+ k ) [y-o  Cx-%-0 _/v 


It  is  no  longer  possible  to  reduce  equation  (4.9)  to  the  form  with  b * 0 by 
changing  the  characteristic  time  as  was  done  in  previous  sections.  The  effect  of 
a phase  lead  component  of  the  nozzle  transfer  function  on  the  stability  of  a 

bipropellant  motor  cannot  be  determined  in  a simple  manner  without  further 

* The  author  is  indebted  to  Prof.  L.  Crocco  for  calling  his  attention  to  the  sign 
mistake  in  equation  (10.1)  in  reference  4 vhich  is  used  in  deducing  equation  (4.9). 


approximation.  Observing  the  fact  that  the  dimensionless  slope  K of  the  adiabatic 
temperature  curve  is  usually  a rather  small  fractional  quantity  except  when  the 
liquid  rocket  is  designed  to  operate  at  very  lean  or  very  rich  mixture  ratio. 

Thus,  as  a rough  approximation,  the  term 

4 jc%-  = k b & j>  <*-«-' ■>-/*<:  (*-*-')] 

may  be  neglected  as  compared  to  other  terms  in  equation  (4.9)  without  affecting 
the  qualitative  trend  of  the  effect  of  the  phase  lead  component  of  the  nozzle 
transfer  function.  After  dropping  this  term  from  equation  "(4.9),  it  can  be  seen 
that  the  new  definition  of  the  characteristic  time  <9^  (l+b)  will  reduce  the 
equation  of  mass  balance  to  the  form  by  putting  b ■ 0.  The  two  equations  of  the 
dynamics  of  the  fuel  and  the  oxidizer  feeding  system,  each  of  them  of  the  fora 
of  equation  (2.10),  can  also  be  reduced  with  the  new  dimensionless  parameters  as 
defined  in  equation  (3.1). 

JL  &e. 

With  this  additional  approximation  of  neglecting  6 qualitative 

conclusion  with  regard  to  the  increasing  stabilizing  effect  of  larger  phase  lead 
component  of  the  nozzle  transfer  function,  as  deduced  In  previous  sections  for  the 
monopropellant  case,  holds  good  for  the  bipropellant  case  as  well  provided  that 
the  feeding  systems  are  not  servo -controlled.  The  stabilizing  effect  manifests 
Itself  as  the  increase  of  the  critical  time  lag  for  neutral  oscillation  in  the 
liquid  rockets.  For  the  bipropellant  case,  however,  the  factor  b of  the  phase 
lead  component  will  in  general  not  only  be  a function  of  the  geometry  of  the 
nozzle  and  the  specific  heat  ratio  y of  the  combustion  gases  as  for  the  mono- 
propellant  rocket,  but  may  also  depend  on  the  magnitudes  of  K and  /e  . The 
entropy  variation  of  the  burned  gas  due  to  the  variation  of  the  stagnation 
temperature  can  no  longer  be  neglected  as  was  in  the  monopropellant  case  where  the 
entropy  variation  is  essentially  due  to  different  dissipative  agents.  Furthermore 
the  neglect  of  the  tens  b which  involves  all  the  feeding  system  parameters, 

may  result  in  a small  dependence  of  the  apparent  value  of  b on  feeding  system  constants . 
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V.  PHYSICAL  ASPECTS  AND  DISCUSSION  OF  RESULTS 

The  fact  that  the  existence  of  a phase  lead  component  of  the  nozzle  transfer 
function  is  in  effect  an  increase  of  the  gas  residence  time  need  some  clarification* 

If  we  examine  the  analyses  (8)  and  (9)  for  the  nozzle  admittance  ratio,  we  see 
that  the  phase  lead  component  arises  because  of  the  inertia  of  the  gas.  The 
inertia  of  the  gas  in  the  subsonic  portion  of  the  nozzle  leads  to  an  accumulation 
of  burned  gas  when  the  pressure  in  the  chamber  is  increased  in  addition  to  the 
accumulation  of  gas  in  the  combustion  chamber*  Therefore  this  constitutes  an 
increase  of  the  effective  capacity  of  the  chamber  in  storaging  the  gas  and  a 
corresponding  increase  of  the  characteristic  time  of  the  system.  In  the  dimensionless 
formulation  of  Crocco  (li)  which  is  followed  in  the  present  treatment,  the  characteristic 
time  is  the  gas  residence  time  defined  as  the  ratio  of  the  mass of  the  burned 
gis  in  the  combustion  chamber  and  the  rate  w of  flow  of  gas  out  of  the  chamber  in 
steady  state  operation, 

s A*  ■ £-*-) 

where  Ve  is  the  volume  of  the  combustion  chamber*  is  the  throat  area  of  the 
nozzle  and  c*  is  the  characteristic  exhaust  velocity.  Vlth  At  and  c*  of  a given 
rocket  fixed,  an  increase  in  the  effective  volume  of  the  combustion  chamber 
corresponds  to  an  increase  in  the  gas  residence  time  • It  is  therefore  natural, 
though  it  may  not  be  quite  right,  to  think  that  the  increase  of  the  characteristic 
time  due  to  the  phase  le^d  component  simply  means  that  the  volume  of  the  subsonic 
portion  of  the  nozzle  should  be  added  to  the  volume  of  the  combustion  chamber  in 
determining  the  characteristic  time.  It  should,  however,  be  noted  that  the  pressure 
and  temperature  level  of  the  gas  in  the  nozzle  is  decreasing  continuously  toward,  the 
throat  and  while  the  chamber  pressure  varies,  the  pressure  in  the  nozzle  also  varies 
but  in  a different  manner.  Therefore  the  capability  of  the  volume  ^ of  the  sub- 
sonic portion  of  the  nozzle  in  accumulating  burned  gas  in 
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response  to  an  Increase  in  the  chamber  pressure  is  different  from  the  capability 
of  cm  equal  volume  9 in  the  combustion  chamber . Knowing  the  nozzle  constant  b 
of  the  phase  lead  component,  the  coefficient  £ of  the  effectiveness  of  the 
nozzle  volume  as  compared  to  chamber  volume  can  be  determined  from 

<V  = e,  (2l+^.L),  e. 


b Oq  ■ 


(5.1) 


For  the  case  of  a nozzle  with  linear  steady  state  velocity  distribution  In  the 
subsonic  portion  and  for  lsentropic  oscillations,  the  coefficient  £ Is  obtained 


from  equations  (2.k),  (2.5)  and  (5.1)  as 

V.«-  -5* ;/ 
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rw-iX/r-^L  -il  (5.2) 

t y+i  ) 3.  ■■ 5/  (fr+J  )*  JJ 

which  is  a function  of  contraction  ratio  $ of  the  nozzle  alone.  Sample 
calculations  for  y * 1.2C  with  Ze.  *>  0.20,  0.10,  and  0.05  corresponding  to  Mach 
number  Me  '*■  O.h,  0.3  and  0.2  respectively  at  the  entrance  section,  the  coefficient 
£ is  equal  to  O.93,  O.89,  and  C.87  for  each  of  the  three  cases.  The  variation 
of  £ for  such  a vide  range  of  entrance  Mach  number  variation  is  rather 
Insignificant  and  the  magnitude  of  £ can  be  approximately  taken  as  0.9. 

Croc co*  has  shown  that  the  admittance  ratio  of  a nozzle  with  velocity 
substantially  linear  upstream  of  the  sonic  throat  but  not  linear  in  the  neighborhood 
of  the  entrance,  can  best  be  approximated  by  the  ratio  of  the  nozzle  where  the 


Unpublished  work  of  Professor  L.  Crocco.  An  outline  of  this  work  will  be 
included  in  the  appendlce'  of  the  volume  on  "Cosibustlon  Instability  in  Liquid 
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velocity  distribution  is  linear  with  velocity  gradient  equal  to  that  at  the 

sonic  throat  of  the  given  nozzle.  Under  such  circumstances,  the  magnitude  of  £ 

would  be  expected  to  be  approximately  0.9.  However,  for  nozzles  whose  velocity 

distribution  is  far  from  being  linear,  we  don't  have  any  idea  as  to  the  relation 

of  b with  the  volume  , and  the  value  of  £ may  be  very  much  different  from  0.9* 

It  cannot  be  over -emphasized  here  to  point  out  that  while  the  concept  of 

increasing  available  chamber  volume  leads  to  a simple  and  convenient  interpretation 

of  the  stabilizing  effect  of  the  nozzle  through  the  increase  of  characteristic 

time)  this  concept,  however,  is  not  established  on  sound  basis  but  simply  a 

conjecture  even  for  low  frequency  oscillations.  The  stabilizing  effect  of  the 

nozzle  toward  high  frequency  oscillation  is  known  to  be  based  on  an  entirely 

different  mechanism  which  has  little  to  do  with  the  volume;  and  the  stabilizing 

effect  of  the  nozzle  toward  low  frequency  oscillations  is  not  clearly  connected 

with  the  volume  of  the  converging  portion  of  the  nozzle.  The  complicated  state 

of  affairs  is  simply  superficially  avoided  by  using  the  coefficient  £ which 

cannot  be  determined  without  the  precise  knowledge  of  the  nozzle  constant  b of 

the  phase  lead  component  of  the  nozzle  transfer  function. 

Having  concluded  the  stabilizing  effect  of  the  factor  b that  determines  the 

phase  lead  component  of  the  nozzle  transfer  function  or  the  deviation  of  the  flow 

in  the  nozzle  from  being  quasi-steady,  we  can  deduce  several  qualitative 

conclusions  concerning  the  effect  of  nozzle  geometry  on  the  lov  frequency  combustion 

stability  of  liquid  rockets.  Analytical  expression  for  b is  known  only  for  the 

« 

special  nozzle  with  linear  steady  state  velocity  distribution  in  the  subsonic 

portion  of  the  nozzle.  It  is  believed,  however,  that  the  change  of  the  velocity 

» 

distribution  near  the  entrance  is  not  going  to  change  radically  the  qualitative 

dependence  of  the  magnitude  of  b on  the  two  parameters  and  ^ i» 

determined  by  tbs  Mach  number  of  the  gas  entering  the  nozzle  and  is  therefore  a 

A * 

parmmster  of  the  contraction  ratio  — , the  ratio  of  the  throat  area  to 

i>  ' o 


18 


the  area  of  entrance  section  of  the  nozzle.  U.x  is  defined  only  when  there  is 
a predominating  slope  in  the  profile  of  the  velocity  distribution.  For  a given 
entrance  Mach  msaiber  and  a given  type  of  velocity  distribution,  essentially 

a parameter  representing  the  "effective  length”  of  the  converging  section  of  the 

nozzle . Since  b is  Increased  by  decreasing  £ e end  , it  is  increased  by 

A* 

decreasing  the  contraction  ratio  or  by  Increasing  the  length  of  the 

converging  section  of  the  nozzle.  Consequently,  we  have: 

1.  Low  frequency  stability  of  a liquid  rocket  is  improved  by  decreasing  the 
contraction  ratio  of  a nozzle  while  keeping  the  length  of  the  converging  section 
constant  and  the  profiles  of  velocity  of  similar  shape. 

2.  Low  frequency  stability  of  a liquid  rocket  is  improved  by  Increasing  the 
length  of  the  converging  section  of  the  nozzle  while  keeping  the  contraction  ratio 
of  the  converging  section  unchanged  and  the  profiles  of  velocity  of  similar  shape. 

The  stabilizing  effect  manifests  itself  analytically  as  an  increase  of  the 
critical  time  lag  for  neutral  oscillation  and  can  be  observed  experimentally  by 
the  decrease  of  the  Intensity  of  "Hugging  and  also  by  the  decrease  of  the  frequency 
of  marginally  unatahl e oscillations  in  the  liquid  rockets. 
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